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Представлено описание математического обеспечения, предназначенного для проведения научных и проектных исследований 
процессов в термоэмиссионных преобразователях, электрогенерирующих элементах и многоэлементных сборках. 


При создании компактных автономных источ- 
ников электрической энергии, обладающих значи- 
тельным ресурсом работы, одним из перспективных 
является термоэмиссионный метод преобразования 
энергии. Данный метод позволяет преобразовывать 
тепловую энергию, полученную с помощью ядерно- 
го источника, непосредственно в электрическую. В 
термоэмиссионных установках протекают многочи- 
сленные взаимосвязанные процессы [1-4]: ней- 
тронно-физические, эмиссионные, плазменные, 
электрические, гидравлические, тепловые, адсорб- 
ционные и др. Экспериментальные исследования и 
испытания термоэмиссионных установок являются 
сложными и дорогостоящими, поэтому значитель- 
ное место в их исследованиях и проектировании за- 
нимает математическое моделирование [3]. 

При разработке математических моделей одной 
из проблем является выбор их типов - стационар- 
ные, нестационарные, статические, эмпирические 
и т.д. Решается эта проблема в зависимости от це- 
лей моделирования на предварительном этапе об- 
следования объекта путем оценки характерных 
времен протекания процессов. Как правило, требу- 
ется провести расчеты статических и динамических 
вольт-амперных характеристик (ВАХ) многоэл- 
ементных термоэмиссионных электрогенерирую- 
щих сборок (ЭГС) при заданных законах измене- 
ния уровней тепловой мощности, давления пара 
цезия, входной температуры и расхода теплоноси- 
теля, а также сопротивления нагрузки. То есть, мо- 
делирование характеристик ЭГС является актуаль- 
ной и сложной задачей, особенно в дуговом режи- 
ме, когда требуется найти взаимосогласованное ре- 
шение нескольких краевых задач, так как при от- 
клонении от равновесного состояния за счет изме- 
нения электронного охлаждения эмиттеров, по- 
тенциального поля ЭГС и джоулева тепловыделе- 
ния в элементах конструкции изменяется темпера- 
турное поле электродов. При расчетах отдельных 
ЭГС обычно влиянием температур на нейтронно- 
физические процессы, а, следовательно, и уровень 
тепловыделения, пренебрегается. Характерное 
время гидравлических процессов зависит от пропу- 
скной способности межэлектродных зазоров 
(МЭЗ), длины ЭГС и изменяется в широких преде- 
лах. Обычно, за исключением проведения спе- 
циальных экспериментов [3], съем ВАХ произво- 
дится после установления давления пара цезия в 
межэлектродных зазорах (МЭЗ), поэтому их инер- 
ционность можно не учитывать. 


Характерное время тепловых процессов соста- 
вляет Гт«10...40 с и определяется теплофизическими 
свойствами используемых материалов, уровнем ге- 
нерируемого тока, условиями охлаждения и размера- 
ми конструкционных элементов установок. В термо- 
эмиссионных установках, в зависимости от состоя- 
ния МЭЗ, может реализовываться вакуумный или 
плазменные (квазивакуумный, прямопролетный, 
диффузионный, дуговой) режимы. Для осуществле- 
ния плазменных режимов в МЭЗ вводят атомы це- 
зия, которые за счет ионизации образуют плазму, а за 
счет адсорбции существенно увеличивают эмисси- 
онную способность электродов. Эмиссионные ха- 
рактеристики электродов ЭГС определяются кине- 
тикой адсорбции (десорбции) атомов цезия. Их ин- 
тенсивность экспоненциально зависит от теплоты 
испарения атомов цезия с электродов. Характерное 
время эмиссии, определяемое временем адсорбции- 
десорбции атомов цезия на электродах, полученное 
на основе оценок [4], составляет г/=1(Г 5 ...1(Г 1 с. Для 
плазменных процессов в термоэмиссионных преоб- 
разователях энергии (ТЭП) существует иерархия ха- 
рактерных времен, которая связана с различными 
подсистемами плазмы. В рабочих условиях оно 
определяется временем релаксации возмущений 
плотности плазмы и составляет т^Ю 6 ... 10 3 с [5]. 
Время установления электрических характеристик 
гораздо меньше, чем для плотности плазмы т э «т п . 

Таким образом, в области времен 5... 50 с тепло- 
вые процессы можно рассматривать в нестацио- 
нарном, а остальные - в квазистационарном при- 
ближении. 

Термоэмиссионные установки можно разделить 
на три группы. В первую группу входят ТЭП - это 
установки, предназначенные в основном для лабо- 
раторных исследований, у которых внешние пара- 
метры (температура эмиттера Т Е , температура кол- 
лектора Т с , давление насыщенных паров цезия Р а , 
межэлектродный зазор сі) практически не меняют- 
ся. Во вторую группу входят электрогенерирующие 
элементы (ЭГЭ) и электрогенерирующие сборки, 
выходные переменные которых существенно зави- 
сят от пространственных координат. На основе ЭГЭ 
и ЭГС проектируются термоэмиссионные ядерные 
энергетические установки третьей группы электри- 
ческой мощностью от 10 до 500 кВт. 

Рассмотрим ЭГС, состоящую из п последова- 
тельно соединенных элементов, упрощенная схе- 
ма, которой представлена на рисунке. 
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Рисунок. Упрощенная схема многоэлементной ЭГС: 1) эмит- 
тер, 2) торцы, 3) ядерное горючее, 4) цилиндриче- 
ская перемычка, 5) дисковая перемычка, 6) дистан- 
ционаторы, 7) коллектор, 8) коллекторный пакет, 
9) токовыводы 


Эмиттер (1) и торцы (2) ЭГЭ представляют со- 
бой оболочку из тугоплавкого материала, как пра- 
вило, вольфрама, внутри которой находится ядер- 
ное горючее (3). Один из торцов эмиттера свобо- 
ден. Ко второму торцу присоединена коммута- 
ционная перемычка, которая состоит из двух ча- 
стей: цилиндрической (4) и дисковой (5). Пере- 
мычка выполнена также из тугоплавкого материала 
- ниобия. Для поддержания заданной величины 
МЭЗ используются дистанционаторы (6), которые 
располагаются по окружности, и их количество мо- 
жет быть различным. Коллектор (7) ЭГЭ предста- 
вляет собой цилиндрическую оболочку из ниобия. 
Сброс тепла производится через коллекторный па- 
кет (8), который в общем случае выполняется в ви- 
де многослойной оболочки из материалов с раз- 
личными теплофизическими свойствами. 

В связи с неравномерностью выделения тепла 
по высоте активной зоны реактора, для поддержа- 
ния одинакового уровня максимальных температур 
эмиттеров в рабочей области используется геоме- 
трическое профилирование ЭГС, которое приво- 
дит к различным длинам ЭГЭ. 

К внешним переменным, определяющим вид 
ВАХ ЭГС, относятся: 

• уровень (ду, Вт/см’) и профиль объемного те- 
пловыделения в сердечниках ЭГС; 

• входная температура теплоносителя ( Т т , К); 

• объемный расход теплоносителя (О с , г/(см 2 с)); 

• давление насыщенных паров цезия (Р с „ Па) в МЭЗ; 

• сопротивление электрической нагрузки ( К н , Ом). 
Зависимость температуры атомов цезия Т с% от 

Р а в резервуаре определяется формулой [4] 


Т = - 


3740 


6 , 78 - 1 ^) 


( 1 ) 


При работе ЭГС тепловая энергия, выделяемая 
в сердечниках, приводит к разогреву эмиттеров и 
торцевых оболочек. Часть энергии преобразуется в 
электрическую за счет эмиссии электронов с эмит- 
тера, а часть сбрасывается тепловым излучением и 
теплопроводностью паров цезия через МЭЗ и кол- 
лектор в систему охлаждения реактора. Через ди- 
станционаторы и коммутационные элементы про- 


исходят потери тепловой энергии с эмиттеров, ко- 
торые приводят к неравномерности температурно- 
го поля электродов. Кроме того, через дистанцио- 
наторы и изолирующие элементы неизбежны утеч- 
ки электрического тока, из-за которых снижается 
величина отдаваемой в нагрузку мощности. Вслед- 
ствие больших величин тока и конечности омиче- 
ского сопротивления эмиттера и коллектора их по- 
верхности не эквипотенциальны. 

В дальнейшем будем использовать цилиндриче- 
скую систему координат с осью і, направленной 
вдоль оси ЭГС слева направо. Начало координат рас- 
положим в плоскости соединения эмиттера первого 
ЭГЭ с токовыводом. Рассматривается азимутально- 
симметричная ЭГС. Под точкой с координатой г по- 
нимается любая точка, расположенная на окружно- 
сти, образованной пересечением эмитирующей по- 
верхности плоскостью, перпендикулярной к оси ци- 
линдра в точке і. В общем случае поверхности эмит- 
тера и коллектора не изотермические и не эквипо- 
тенциальные, поэтому работы выхода электродов, 
разность потенциалов между электродами, потери 
тепловой энергии с эмиттера, следовательно, и плот- 
ность генерируемого тока эмитирующей поверхно- 
сти являются функциями координаты і. Кроме того, 
вывод газообразных продуктов деления топлива в 
МЭЗ ухудшает эмиссионные свойства электродов и 
повышает теплопроводность зазора. 

При прохождении электрического тока по 
электродам в зазоре образуется собственное маг- 
нитное поле, следствием которого является возни- 
кновение двух эффектов [1,6]: 

• уменьшение подвижности электронов, а, следо- 
вательно, и тока насыщения преобразователя в 
диффузионном режиме; 

• появление пондеромоторных сил, которые соз- 
дают градиент давления в закрытых, либо дви- 
жение плазмы в открытых ЭГС. 

В течение ресурсных испытаний к перечислен- 
ным особенностям объекта исследования следует 
отнести деформацию топливно-эмиттерных бло- 
ков, а также изменения эмиссионных, тепло- и 
электрофизических свойств материалов. 

То., плотность генерируемого тока оказывается 
сложной функцией многих переменных в каждой 
точке координаты г, а именно, температур электро- 
дов, давления паров цезия, напряжения между 
электродами, межэлектродного зазора и т.д. 

Для получения этой зависимости могут быть 
использованы экспериментальные данные, эмпи- 
рические [7, 8] или теоретические [1, 9] модели 
ВАХ ТЭП. В основу модели, изложенной в [7], по- 
ложены обнаруженные закономерности зависимо- 
сти внутреннего падения напряжения на зазоре от 
параметров преобразователя, которые выражаются 
следующим образом: 


К = ѵ с1 (Р Сі сі,т Е ,.;)+Аѵ а 


кТ 

р Сі і,— іп(./ Ё Ц) 

е 


, ( 2 ) 
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где к - постоянная Больцмана; е - заряд электро- 
на; ^ - ток, проходящий через преобразователь; 
- эмиссионный ток эмиттера. 

Из (2) видно, что асимптотическая составляю- 
щая падения напряжения на зазоре Ѵ іж зависит от 
Р,/І, Т Е и /. Асимптотическая ветвь ВАХ реализует- 
ся в дуговом режиме работы при условии К<І Е . В 
области токов /— >/ г падение напряжения на зазоре 
значительно возрастает за счет дополнительной со- 
ставляющейДИ Л зависящей от Р а Ви (кТ І /ё)\а.Ц 1 /І). 

Выходное напряжение преобразователя рассчи- 
тывается по формуле: 

Г н ={Р?-Р с )!е-Г,. ( 3 ) 


Эффективная работа выхода эмиттера Р Е Ф берет- 
ся как наилучшая оценка из трех величин {Р в ,Р Е ,/и р } 


Р и < р < р 

1 Е’ Ир ^ 1 Е 


^ ф = 


Рр, Р Е<Ир< Р В 
^В’^Е < ^В < Рр 


Р Р > Р 

1 В’ 1 Е ^ * В 


(4) 


где Р в ,Р Е ,ц р - «истинная» работа выхода, больцманов- 
ская работа выхода и химический потенциал плазмы, 

Р в =кТ Е \п(тТ Е 2 /Р), (5) 


Р р =1,945 +1,15*7И2,5ВД) -1§(^|) -3,0]. (6) 

Работы выхода электродов определяются с по- 
мощью кривых Рейзора [4] Рцо Т Е (о / Т Сі , Рщо), 

где Р° т - вакуумная работа выхода эмиттера (кол- 
лектора), Т 0 . определяется по формуле (1). Эмисси- 
онные плотности тока электродов определяются по 
формуле Ричардсона - Дэшмана [4] 

Л(С) = 120 Те(С) еХ Р( (С ) /( ^ (С ) )) • 

В [7] также показано эффективное применение 
модели для расчета ВАХ ТЭП, не вошедших в се- 
мейство экспериментальных характеристик, ис- 
пользованных при ее построении. Эмпирическая 
модель позволяет проводить расчеты ВАХ при из- 
менении параметров ТЭП в следующих диапазонах: 

7/=1600...2100 К, 6</<30 А/см 2 , 
й^ 0,01 см, 0,4788<Р & /<83,79 Паюм. 

При построении модели использовались экспе- 
риментальные ВАХ ТЭП с пренебрежимо малой 
эмиссией с коллектора. Составляющая падения на- 
пряжения на зазоре ДКД_Р с Д(&7//е)1п(/ г /./)], кото- 
рая учитывает эффект Шоттки, определена только 
для ///>1. Так как при существенно возра- 

стает эффект Шоттки, то для учета «горячего» кол- 
лектора в рассматриваемую модель внесены сле- 
дующие усовершенствования. 

Во-первых, т.к. значения тока, при которых про- 
изводится расчет напряжения между электродами, 
ограничиваются областью до точки перегиба ВАХ 
/<^ пер (Р Сі /)4 то в области более высоких токов ха- 
рактеристика аппроксимируется прямой линией, 
параметры которой так же рассчитываются на ос- 


нове эмпирической модели в точке перегиба. Ко- 
эффициент к пср определяется по формуле, аппрок- 
симирующей экспериментальные данные [7]. 

Во-вторых, влияние электронной эмиссии с 
коллектора учитывается за счет дополнительного 
падения напряжения на зазоре 

АѴ =1°’ ССЛИ ' У ' : ~ кж ’ ] ' 

с \кТ с / е 1п(У с /(& ое У)), если / с > & ое /, (7) 

где / с , к ж - эмиссионный ток с коллектора и коэф- 
фициент, учитывающий эмиссию с коллектора. Ре- 
комендуемое значение к о =0,2. 

В-третьих, область диффузионного режима ра- 
боты ТЭП аппроксимируется зависимостью вида 

Ѵ н (У) = Ѵ# — кТ Е 1п( У* / У ) / е, (8) 

где Ѵ Е , Р - минимальные напряжение и ток дуго- 
вого режима работы ТЭП. Эти значения определя- 
ются как координаты точки пересечения дуговой 
ветви ВАХ ТЭП и горизонтальной прямой, отсе- 
кающей по оси / значение, равное току насыщения 
в диффузионном режиме [ 1 ] . 

Описанные модификации модели, уравнения 
(2-8), позволили существенно расширить область 
расчета ВАХ ТЭП, что очень важно при проведе- 
нии многовариантных расчетов характеристик 
ЭГС в автоматическом режиме. 

С помощью модифицированной модели [7] 
можно рассчитать набор ВАХ ТЭП в плазменных 
режимах. В [9] приведены модели, позволяющие 
рассчитывать характеристики преобразователя в 
вакуумном и квазивауумном режимах. 

То., для исследования характеристик ЭГС зави- 
симость плотности генерируемого тока от параме- 
тров преобразователя можно получить тремя спо- 
собами. Во-первых, использовать эксперименталь- 
ные характеристики в виде таблиц. Во-вторых, 
включить модель ВАХ ТЭП в модель расчета ВАХ 
ЭГС. В-третьих, по модели расчета ВАХ ТЭП рас- 
считать зависимости плотности генерируемого то- 
ка от внешних параметров в виде таблиц и далее их 
использовать при моделировании ЭГС. Недостат- 
ком второго способа является необходимость мно- 
гократного расчета одних и тех же характеристик 
ТЭП, что приводит к большим затратам времени. 
Поэтому на практике наибольшее применение на- 
шли первый и третий способы. 

Для расчета статических ВАХ коаксиальных 
ЭГС к настоящему времени разработано большое 
количество моделей [10-14]. Все они базируются 
на решении с помощью ПЭВМ подобных систем 
дифференциальных уравнений, описывающих те- 
пловые и электрические процессы в ЭГЭ, а разли- 
чаются методами расчета и разной степенью дета- 
лизации учета особенностей их конструкций. Эти 
различия определяют выбор математической моде- 
ли для решения конкретной задачи. Так для оце- 
ночных расчетов вполне пригодны простые моде- 
ли, позволяющие оперативно рассчитывать ВАХ 
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ЭГС [10]. В то же время при изучении процессов в 
ЭГЭ и ЭГС, а также при прогнозировании ожидае- 
мых характеристик и их применении для решения 
задач идентификации целесообразно использова- 
ние детальных математических моделей ЭГС 
[11-14]. Под детальностью модели понимается 
максимально возможная степень отражения схе- 
мно-конструкционных и технологических особен- 
ностей ЭГС и изменяющихся условий работы. 

Следует отметить, что одним из наиболее ин- 
формативных режимов работы ЭГС, позволяющим 
получить дополнительную информацию о ее со- 
стоянии, является динамический. Так на основе 
обработки экспериментальных переходных про- 
цессов разработан ряд методов оценки внутренних 
неизмеряемых параметров, а с помощью модели 
[12] объяснено аномальное поведение некоторых 
ЭГС [3]. 

Так как тепловые процессы более инерционны 
по сравнению с электрическими, то математиче- 
скую модель ЭГС можно представить в виде сово- 
купности нестационарной тепловой и стационар- 
ной электрической моделей [12]. 

При составлении математической модели те- 
пловых процессов ЭГС были приняты допущения 
кроме азимутальной симметрии, о пренебрежении 
перепадом температур по толщине электродов, 
коммутационных перемычек и несущей трубки. 

Уравнения, описывающие тепловые процессы в 
элементах конструкции ЭГС, записываются следу- 
ющим образом: 


д_ 

дг 


а,(т,) 


дТ і 

дг 


)+' д (гА 1 (Т) д %+ Чѵі ( г ,г) = 
г дг дг 

= с І (Т І )р,(тА 

ді 


(9) 


где Я;, с„ рі - теплопроводность, удельные теплоем- 
кость и плотность материала /-ого элемента кон- 
струкции; д ѴІ ( 7 „г) - удельное тепловыделение в /-ом 
элементе конструкции. 

Уравнение теплопроводности жидкометалличе- 
ского теплоносителя: 


А (Я г (Т т ) д ^.) + С с с т {Т т ) д ^+ <* т (г) = 

02 02 02 

= с т {Т т )р т (Т т )^~. ( 10 ) 
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тепловой поток с эмиттера 

-Я, (Т Е ) ^ = ](сі, Р а ,Т Е ,Т с ,Щг))х 
дг 


х[Ѵ{2) + Чс (Т с ,Т С5 )Іе\ + 


+е (Т Е , Т с )ст (Т Е -Т с 4 )+ Х ^ Г Е’ Т с’ Р а) 


( 12 ) 


ІТ Е -Т С ), (13) 


где Т г =Т(г=г г ), Т с =Т(г=г с ), г Е , г с - наружный ради- 
ус эмиттера и внутренний радиус коллектора соот- 
ветственно, в - приведенная степень черноты из- 
лучающих поверхностей, ст - постоянная Стефана- 
Больцмана, < 7 С - теплота конденсации электронов 
на коллекторе; 

- тепловой поток на коллектор 
дТ 

-Я с (Т с ) с = )Ы, Р а ,Т Е ,Т С , Ѵ(2)) X 
дг 


*[<7 С (Т с ,Т а )/ё\ + е(Т Е , Т с )а(Т 4 - Т с 4 ) + 

+ я с ,(т Е ,т с ,р с ,) { т е - Тс) . 


сі 


(14) 


тепловой поток между поверхностями, распо- 
ложенными на расстоянии 8і, 


ЦтА го =е(Т)а(Т 4 

дг 


-Т 

20 ± 20+5г 


) + 


*ъ(Т) 


(Т 20 Т 20+52 )> 


(15) 


- тепловой поток на несущую трубу 

-1 с (Т) дТс =а(Т)(Т с -Т ІІГ ), 
дг 

где а - коэффициент теплопередачи коллекторно- 
го пакета; 


(16) 


входная температура теплоносителя 

Т т \ы=Т Ій . (17) 

Начальные условия задаются следующим образом 


ТО, гр) =Т(2,Г,1 0 ), 


(18) 


Уравнения (9, 10) дополняются граничными 
условиями, которые имеют следующий физиче- 
ский смысл: 


- равенство нулю теплового потока на оси в силу 
радиальной симметрии ЭГС 


дТ 

дг 



(И) 


- равенство температур и тепловых потоков на 
границах конструкционных материалов 


где Г(?,г,/ 0 ) - температурное поле ЭГС в стацио- 
нарном состоянии. 

Распределение потенциалов по длине электро- 
дов к-ото ЭГЭ сборки с учетом плотностей токов 
утечки между электродами и через коллекторную 
изоляцию ] ш описывается уравнениями вида [12] 
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Технические науки 


На границах к-ото ЭГЭ с перемычками должны 
выполняться следующие условия 
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где 5” - площади поперечных сечений соответ- 
ствующих элементов ЭГЭ; р - удельные электри- 
ческие сопротивления материалов; Т Р - температу- 
ра перемычки; Г ь Г 2 и Г 3 - границы соприкоснове- 
ния электродов с коммутационной перемычкой. 

Напряжение на нагрузке находится по формуле 


Ѵн = (У с 


2 = 0 


-К 


) 4і^Пі (23) 


причем, ІІ н =І к К н К т , К п - сопротивления токовы- 
водов, І к , І ш - токи, проходящие по соответствую- 
щим токовыводам. 

Для расчета динамических характеристик ЭГС 
необходимо найти самосогласованное решение 
уравнений (9, 10, 19) с граничными (11-17, 20-23) 
и начальными (18) условиями при внесении раз- 
личных возмущений в математическую модель 
сборки. 

Таким образом, представленное математиче- 
ское описание процессов, происходящих в ТЭП и 
ЭГС, дает возможность решить ряд текущих задач, 
возникающих как на этапе проектирования, так и 
на этапе испытаний ЭГС, а именно, рассчитывать: 

• статические и динамические ВАХ ЭГС; 

• переходные процессы; 

• аварийные ситуации (короткое замыкание, 
обрывы перемычек, ухудшение свойств коллек- 
торной изоляции и т.п.); 

• длины ЭГЭ в ЭГС при геометрическом профи- 
лировании; 

• влияние различных параметров на выходные 
характеристики ЭГС; 

• внутренние параметры ЭГС по выходным ста- 
тическим и динамическим характеристикам и 
т.д. 
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